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Abstract 

Adaptation of the chemotaxis sensory pathway of the bacterium Escherichia coli is integral for detecting 
chemicals over a wide range of background concentrations, ultimately allowing cells to swim towards 
sources of attractant and away from repellents. Its biochemical mechanism based on methylation and 
demethylation of chcmoreccptors has long been known. Despite the importance of adaptation for cell 
memory and behavior, the dynamics of adaptation arc difficult to reconcile with current models of pre- 
cise adaptation. Here, we follow time courses of signaling in response to concentration step changes 
of attractant using in vivo fluorescence resonance energy transfer measurements. Specifically, we use a 
condensed representation of adaptation time courses for efficient evaluation of different adaptation mod- 
els. To quantitatively explain the data, we finally develop a dynamic model for signaling and adaptation 
based on the attractant flow in the experiment, signaling by cooperative receptor complexes, and multiple 
layers of feedback regulation for adaptation. We experimentally confirm the predicted effects of changing 
the enzyme-expression level and bypassing the negative feedback for demethylation. Our data analy- 
sis suggests significant imprecision in adaptation for large additions. Furthermore, our model predicts 
highly regulated, ultrafast adaptation in response to removal of attractant, which may be useful for fast 
reorientation of the cell and noise reduction in adaptation. 

Author Summary 

Bacterial chemotaxis is a paradigm for sensory systems, and thus has attracted immense interest from bi- 
ologists and modelers alike. Using this pathway, cells can sense chemical molecules in their environment, 
and bias their movement towards nutrients and away from toxins. To avoid over- or understimulation of 
the signaling pathway, receptors adapt to current external conditions by covalent receptor modification, 
ultimately allowing cells to chemotax over a wide range of background concentrations. While the ro- 
bustness and precision in adaptation was previously explained, we quantify the dynamics of adaptation, 
important for cell memory and behavior, as well as noise filtering in the pathway. Specifically, we study 
the intracellular signaling response and subsequent adaptation to concentration step changes in attrac- 
tant chemicals. We combine measurements of signaling in living cells with a dynamic model for strongly 
coupled receptors, even including the effects of concentration flow in the experiment. Using a novel way 
of summarizing time-dependent data, we derive a new adaptation model, predicting additional layers of 
feedback regulation. As a consequence, adaptation to sudden exposure of unfavorable conditions is very 
fast, which may be useful for a quick reorientation and escape of the cell. 
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Introduction 

Cells are able to sense and respond to various external stimuli. To extend the working range of their 
sensory pathways, biochemical mechanisms allow for adaptation to persistent stimulation, resulting in 
only a transient response. The dynamics of adaptation are important as they often represent the cells' 
memory of previous environmental conditions, directly affecting cellular behavior [IHZ]- Fast adaptation 
means that cells stop responding and that their biochemical pathways quickly "forget" the stimulus. In 
contrast, slow adaptation leads to long-lasting effects in the cells. The dynamics of adaptation are often 
difficult to understand in detail, since they emerge from multiple, simultaneously occurring processes. 
In higher organisms, adaptation is best documented in the insect and vertebrate visual system, where 
multiple mechanisms adjust the receptor sensitivity to ambient light levels. For instance, phototransduc- 
tion in the vertebrate eye involves up to nine different mechanisms for adaptation [Q. However, even 
in the well-characterized chemotaxis sensory system in Escherichia coli [51-113] . adaptation, in particu- 
lar its molecular mechanism and dynamics, is not well understood. This constitutes a huge deficit as 
there has recently been immense interest in the chemotactic behavior of bacteria and noise fil- 

tering [TTlQjlJUn] • Here, we use adaptation time-course data from in vivo fluorescence resonance energy 
transfer (FRET) measurements and quantitative modeling to address this problem. 

The chemotaxis pathway in E. coli allows cells to sense chemicals and to swim towards more favorable 
environments by successive periods of straight swimming (running) and random reorientation (tumbling). 
Transmembrane chemoreceptors, including the highly abundant Tar and Tsr receptors, cluster at the cell 
poles and act as "antennas" to detect various attractants (e.g. certain amino acids and sugars) and 
repellents (e.g. certain metal ions) with high sensitivity |21) . Receptors activate an intracellular signaling 
pathway, which results in the phosphorylation of diffusible response regulator CheY (CheY-P) via kinase 
ChcA. CheY-P binds to the flagellated rotary motors and induces tumbling. For details of the pathway 
see the Supplementary Text SI. The interactions between different proteins in the chemotaxis pathway 
during signaling have been well characterized. Specifically, FRET measurements on the response regulator 
CheY-P and its phosphatase CheZ have elucidated the signaling in the chemotaxis pathway [2"2"H2~i] . 

Adaptation in E. coli is based on reversible methylation and demethylation of receptors at specific 
modification sites, catalyzed by enzymes CheR and phosphorylated CheB (CheB-P), respectively. Tar 
and Tsr receptors have four major modification sites. In addition, the Tsr receptor has two minor modifi- 
cation sites which are methylated less strongly [25] . Receptor modification regulates the receptor activity 
and provides a recording of experienced concentration changes 10,26,27 . As a consequence, the rate of 
tumbling was found to adapt precisely for different ligand concentrations [281129] . To achieve the return of 
the receptor activity to its pre-stimulus value, receptor activity-dependent phosphorylation of CheB pro- 
vides a negative feedback on the receptor activity. In addition, the rates of methylation and demethylation 
depend on the receptor activity [30l432j . representing further layers of feedback regulation. To modify 
receptors, CheR and CheB molecules can bind to a specific tether sequence at the carboxyl-terminus of 
Tar and Tsr receptors, and act on several nearby receptors, so-called assistance neighborhoods [33]. This 
is believed to compensate for the low numbers of CheR and CheB (hundreds of molecules) [34] , although 
larger numbers have been reported [35] . 

Although a lot is known about the components of the chemotaxis pathway, several open questions 
remain to be answered in adaptation, (i) Despite their importance for cell behavior, memory and noise 
filtering, the dynamics of adaptation and the methylation level are largely unknown. This is because 
the methylation level is difficult to measure precisely, relying on quantification of receptor protein and 
radioactively-labeled methylation substrate (methionine) incorporated into receptors [25 ]l36ti38] , So far, 
only the initial rate of adaptation was inferred from the rate of change in motor bias in response to 
saturating amounts of added attractant [29] . (u) The molecular mechanism of adaptation, in particular 
demethylation, remains unclear. While CheR binds strongly to the tether, suggested to increase its con- 
centration in the vicinity of methyl-accepting sites [39] , the binding affinity of CheB was found to be very 
low [40]. Instead, binding of ChcB-P to the tether induces an allosteric activation of the receptor, increas- 
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ing the demethylation rate [40j . Furthermore, while the receptor activity-dependence of the methylation 
and demethylation rates is believed to be a requirement for robust precise adaptation (see below) , it is not 
known if adaptation is precise at the receptor level. Time-course data from in vivo FRET experiments, 
monitoring receptor activity upon stimulation, is ideally suited to study the adaptation dynamics and 
address these questions. 

Extensive mathematical modeling has described different aspects of the chemotaxis pathway. However, 
modeling has mainly focused on explaining the initial response to addition of attractant, as well as precise 
adaptation, i.e. the complete return of the signaling activity to prc-stimulus level long after the stimulus. 
Briefly, the Monod-Wyman-Changeux (MWC) model was used to successfully describe the signaling 
of two-state receptor complexes, formed by 10-20 strongly interacting receptor dimers [241I4TH44] . In 
this model, receptor-receptor coupling provides a mechanism for signal amplification and integration. 
Alternative receptor models are outlined in the Discussion. Furthermore, Barkai and Leibler showed 
that precise adaptation is robust (insensitive to parameter variations in the pathway), if the kinetics 
of receptor methylation depends only on the activity of receptors and not explicitly on the receptor 
methylation level or external chemical concentration [45| . Their idea was later extended by others, 
providing conditions for precision [46l|47], as well as robustness to noise by the network architecture [48] 
and assistance neighborhoods ,42, 49]. Most recently, adaptation to exponential ramps and sinusoidal 
concentration changes was investigated [20] . However, none of these studies have directly compared to 
adaptation time-courses from FRET. 

Here, we use in vivo FRET data obtained from cells adapted to ambient concentrations of attractant 
a-methylaspartate (MeAsp; a non-metabolizable variant of amino acid aspartate) and stimulated in a flow 
chamber by various concentration step changes [23]. Recording the average initial response amplitudes 
for each added and, after adaptation, removed concentration step change results in dose-response curves 
(Fig. [l] symbols). We use a dynamic version of the MWC model, which, in addition to mixed complexes 
of Tar and Tsr receptors, includes a more detailed description of the adaptation dynamics than used 
in previous models of chemotaxis. Specifically, we predict multiple layers of feedback regulation during 
adaptation, especially for demethylation by CheB. In addition, we take into account the kinetics of 
attractant flow in FRET experiments. This allows us to quantitatively describe dose-response curves 
(Fig. [T] lines), in particular the observed reduced response amplitudes for removal of MeAsp, which 
previously could not be explained by the MWC model {Inset). To analyze the adaptation dynamics, 
we use the data collapse, a condensed representation of time courses. This data collapse enables us 
to evaluate the effect of ligand flow and adaptation imprecision, to infer the kinetics of the receptor 
methylation level, as well as to efficiently compare adaptation models from the literature to experimental 
data. Finally, we experimentally test two predictions to validate our adaptation model. We change the 
adapted receptor activity, and use a non-regulatable CheB mutant to bypass its negative feedback on the 
receptor activity. Our combined study of experiments and modeling shows that adaptation is relatively 
imprecise at the receptor level for large stimuli, and that demethylation is more tightly regulated than 
previously thought. This leads to very short tumbles for sudden occurrences of unfavorable conditions, 
allowing cells to quickly reorient their swimming direction after a short tumble. 

Results 

Dynamic MWC model for in vivo FRET response 

Our dynamic MWC model, described in the following, combines the previously used MWC model for 
receptor signaling by strongly-coupled receptor complexes (denoted here by static MWC model), with the 
dynamic effects of adaptation by receptor modification, as well as ligand concentration flow. In the static 
MWC model, mixed receptor complexes composed of Tar (aspartate receptor) and Tsr (serine receptor, 
which also binds aspartate with lower affinity) are considered in their in vivo ratio. Using a two-state 
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assumption, the activity of a receptor complex is given by its probability to be in on (active), which 
depends on the free-energy difference F = F on — F Q g between its on and off (inactive) state [HIES : 

e -F on ^ 1 

A = p on = 5 ; p — = - — ; p. (f) 

This free-energy difference, F(m,c), is determined by two contributions, one from methylation (in terms 
of receptor methylation level m) favoring the on state, and one- from attractant binding at MeAsp con- 
centration c favoring the off state. The free-energy difference F also depends on several parameters such 
as free-energy difference per added methyl group, the number N of receptor dimers in a complex, as well 
as the ligand dissociation constants K°f s s and for Tar (Tsr) receptors in their on and off states, 

respectively Most of these parameters were determined previously (see Materials and Methods). Similar 
free-energy based two-state models were recently used to describe clustering of ion channels [50| and small 
GTPases in eukaryotic cells [ST]. In the new dynamic MWC model, we include the effects of variable 
receptor complex sizes, adaptation dynamics, and MeAsp concentration flow on the initial response to 
concentration changes. 

The dependence of the receptor complex size on the ambient concentration and hence methylation 
level was determined as follows: First, the receptor complex size was obtained for each ambient con- 
centration using a least-squares fit to addition dose-response curves (see Fig. [2]4. and Materials and 
Methods). Consistent with previous modeling results, we find that the receptor complex size increases 
with increasing ambient concentration [411152] . As the simplest assumption, we used a linear relationship 
between receptor complex size and ambient concentration (Fig. HJfV), allowing us to interpolate the recep- 
tor complex size for removal dose-response curves. Analyzing the signaling pathway of E. coli, we also 
found the phosphorylation reactions are sufficiently fast to assume that concentrations of phosphorylated 
(and unphosphorylatcd) proteins are in quasi-steady state. Furthermore, the concentrations of activated 
proteins are approximately proportional to the receptor complex activity. Both these conditions allow 
us to use the receptor complex activity as a substitute for the down-stream activity measured by FRET 
reducing the number of model parameters for fitting to data greatly (see Supplementary Text SI). This 
approximation was also used in previous work, but was never explicitly tested [411443] . 

Adaptation occurs on a similar time scale as the kinetics of the MeAsp concentration flow. In ex- 
periments, changes in MeAsp concentration are established over several seconds, due to the finite flow 
speed and mixing effects in the flow chamber. In our model, we assume exponentially rising and falling 
concentration changes upon addition and removal in line with previous measurements by Sourjik and 
Berg (Fig. [2J3) [23]. Adaptation is mediated by methylation and dcmcthylation enzymes CheR and 
ChcB, respectively. The process is described by the kinetics of the average receptor methylation level m 
in a receptor complex, 

dm , .. ,i 

— = 9R (1-A)- 9B A 3 , (2) 

where the adapted receptor-complex activity A* is determined by the steady-state condition dm/dt = 
= <?_r(1 — A*) — gsA* 3 . According to our model, receptors are methylated when the complex is 
inactive, and demethylated when it is active. Furthermore, we postulate a very strong sensitivity of 
the demethylation rate on activity, possibly due to cooperativity of CheB-P molecules. This mechanism 
explains the strong asymmetry, which is observed in experimentally measured time courses (cf. Fig. |2p) 
where adaptation of inactive receptors (methylation) is slow compared to the rapid adaptation of active 
receptors (demethylation). Hence, during a concentration step change the initial response amplitude of 
receptor complexes is reduced by simultaneous adaptation, which is particularly important for removal 
of concentration (see Fig. [2J3 Inset). Note that the asymmetry between slow adaptation of inactive and 
active receptors, respectively, cannot simply be changed by adjusting the rate constants of methylation 
and demethylation individually, since they are constrained by the adapted activity A* . For details of this 
adaptation model see Materials and Methods, and for a potential molecular mechanism of demethylation, 
sec Discussion. 
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Experimental dose-response curves (Fig. [I] symbols) describe the initial response of adapted wild- 
type cells to sudden changes (addition and removal) in McAsp concentration |23j . These responses are 
taken from time courses measured by in vivo FRET (cf. Fig. [5]). Additional, previously unpublished 
dose-response curves are provided in the Supplementary Text SI. For details of the experiments see 
Material and Methods. Our dynamic MWC model, which includes the effects of adaptation and MeAsp 
flow, quantitatively describes the experimental dose-response curves. Specifically, adaptation leads to 
a non-saturated response for large MeAsp step changes Ac at high ambient concentrations, which is 
not seen in the static MWC model without adaptation dynamics (Fig. [1] Inset). Note, however, that 
responses eventually do saturate according to the dynamic MWC model for even larger concentration 
step changes due to the presence of Tsr receptors (at 0.5 mM ambient for Ac>40mM; not shown). The 
dynamic MWC model describes the dose-response data significantly better than the static MWC model, 
as indicated by their overall squared errors in the caption of Fig. [TJ as well as residual errors detailed 
in the Supplementary Text SI. The receptor-complex activity, as well as FRET data were normalized 
by their adapted prc-stimulus values at ambient concentration to compare model and experimental data 
(see Materials and Methods). 



Data collapse of time courses for adaptation dynamics 

The short-term behavior in the time-course data (Fig. [2jU) is essential in deriving our adaptation model, 
used to quantitatively describe dose-response curves (Fig. [I]). Can our adaptation model also describe the 
long-term behavior in the time-course data, and hence the complete adaptation dynamics? Our model 
for precise adaptation predicts that the observable rate of activity change is given by 

dA _ dAdm + dA dc . 
dt dm dt dc dt ' 

where the rate of change of the methylation level dm/dt is described by Eq. [31 and dc/dt is the rate 
of change of the MeAsp concentration. After a concentration step change, the MeAsp concentration is 
constant with dc/dt=0, and the rate of activity change is given by 

dA dAdm .N r 3 , 

~dt = d~m~~dt = ~ ~2~ ^ " A) ~ ^ S f ^ (4) 

where we used that dA/dm=(dA/dF)(dF/dm)=A(l — A)N/2 (see Material and Methods). Hence, the 
rate of activity change is a function f(A) of the activity only, independent of ligand concentration and 
receptor methylation level (except for the minor dependence of the receptor complex size on the ligand 
concentration, see Supplementary Text SI). This predicts a data collapse of all adaptation time courses, 
independent of the duration, size and number of concentration step changes, onto a single curve dA/dt= 
f(A) (Fig. [3]^, thick gray line). This non-monotonous function of the activity has three fixed points: the 
adapted activity A=A*, where methylation and demethylation rates exactly balance each other, as well 
as A~0 and A=l, where the receptor complex activity is saturated in the o/f and on state, respectively. 
Figure |3]?V Inset shows the experimental rate of activity change as extracted from our quantitative time- 
course data from FRET for different concentration step changes at an ambient concentration. We observe 
that, in contrast to the prediction of the model, the rate of activity change depends on the magnitude of 
the concentration step changes. For addition of large concentration step changes (blue symbols), the rate 
is reduced and the activity stays below the pre-stimulus value. Furthermore, for total removal of MeAsp 
concentration (replacement with buffer medium, green symbols), the magnitude of the rate is reduced 
and the activity remains above the pre-stimulus value. 

To explain the deviations from the predicted data collapse, we consider the effects of MeAsp flow and 
imprecise adaptation in our model. According to Eq.[31 each of the two effects contribute independently to 
the rate of activity change. First, we include the MeAsp flow for concentration step changes as described, 
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and simulate time courses based on the precise adaptation model (Fig. solid lines). We find that in 
the demethylation regime (negative rate of activity change), the kinetics of concentration step removal 
gives rise to minor deviations from the curve f(A) in qualitative agreement with experiment. However, in 
the methylation regime (positive rate of activity change), unlike the experimental data, all time courses 
lie accurately on the f(A) curve. Next, we consider imprecise adaptation, i.e. the incomplete return of the 
activity to pre-stimulus level, which is apparent in the time courses (Fig. [2]C and Supplementary Text SI 
for quantification). In our model for imprecise adaptation, Eq. [7] in Materials and Methods, the kinetics 
of the methylation level dm/dt depends explicitly on the receptor methylation level, which leads to sig- 
nificant deviations from the data collapse (Fig.[5]4., dashed lines). Adaptation after addition of increasing 
concentration step changes results in a reduced adapted receptor complex activity (adapted activity after 
removal is always the same as the concentration is the ambient concentration) . Total removal of MeAsp 
concentration (buffer) results in an increased adapted activity. Our imprecise adaptation model is in line 
with the experimental data, showing that the data collapse is an effective way to compare experimental 
and theoretical time courses and to quantify the effects of ligand flow and imprecise adaptation. We also 
studied the effect of changes in receptor-complex size on the data collapse, which we found to be minor 
for the concentrations considered here (see Supplementary Text SI). 

In addition to the adaptation dynamics, the data collapse allows us to determine the kinetics of the 
receptor methylation level, which is difficult to measure directly. Figure [3)3 shows the rate of change of 
the methylation level as a function of the receptor complex activity for experimental data, as well as the 
dynamic MWC model. The data and curves were obtained by dividing the rate of activity change dA/dt 
following concentration step changes by A(l — A). If the activity change is caused only by the adaptation 
dynamics, this procedure yields a function proportional to the rate of change of the methylation level, 
dm/dt. According to our precise adaptation model Eq. [2j the rate of change of the methylation level 
is a monotonically decreasing function of activity with one steady state, marking the adapted receptor 
complex activity (Fig. [3}3, thick gray line). Corresponding to the rate of activity change in Fig.[3]4., the 
kinetics of ligand flow upon concentration step changes, as well as imprecise adaptation result in deviations 
from this curve. As before, we mainly find signatures of imprecise adaptation in the experimental data 
in Fig. |3E Inset. 

Comparison of different adaptation models 

The data collapse of experimental time courses enables the efficient evaluation of different adaptation 
models, including our model and other models from the literature (Fig.0J\). All models considered here 
show precise adaptation with the rates of methylation and demethylation only depending on the receptor 
complex activity, and the explicit activity dependencies given respectively by the first and second term 
in the legend of Fig. [4] For instance, the first model (solid red line) is given by Eq. [2] Two classes of 
models are analyzed here. The first class of models, including our model, is based on a linear activity- 
dependence of the methylation and demethylation rates for binding of CheR and CheB to inactive and 
active receptor, respectively. Feedback from the activity-dependent phosphorylation of CheB is accounted 
for by additional factors of the receptor complex activity. Our model includes cooperative CheB feedback 
(solid red line), while linear CheB feedback (dashed red line) and no CheB feedback (dotted red line) are 
considered as well [15 ,42,49,53]. Another class of models has been proposed, showing ultrasensitivity with 
respect to CheR and CheB protein levels. In these models, the activity- dependence of the methylation 
and demethylation rates for enzyme binding is described by Michaelis-Mentcn kinetics, and linear CheB 
feedback (solid blue line) and no CheB feedback (dashed blue line) is considered [17]. The last model 
has the property that the rate of change of methylation level becomes independent of activity around the 
steady-state, leading to extremely long adaptation times. Details of the alternative adaptation models 
and the fitting procedure are given in the Supplementary Text SI. While several models are consistent 
with the experimental data, our model compares most favorably. The ultrasensitive Michaelis-Mentcn 
model without CheB feedback seems not to be consistent with the data. Comparing simulated time 
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courses for the different adaptation models in Fig. 0J3, our model is best to capture the experimentally 
observed asymmetry between adaptation to addition and removal of concentration step changes. The 
quality of fit between the respective models and data is indicated by their least-squares errors in the 
caption of Fig. @] 

Predictions 

To further validate our adaptation model, we experimentally tested two predictions. First, in our precise- 
adaptation model the data collapse depends strongly on the steady-state activity. For instance, increasing 
the steady-state activity from A*ml/3 to 1/2 changes the data collapse from the solid to the dashed red 
line in Fig. UK. Such an increase in the steady-state activity can be achieved by decreasing CheB expression 
level, corresponding to a decreasing dcmcthylation rate, at constant ChcR expression level. To validate 
this prediction, a different wild-type strain (WT2) was created, in which CheB expression was induced 
from a plasmid, while all other chemotaxis proteins were expressed as before (WT1). The steady-state 
activity was estimated to be A* ml/ 2 (compared to 1/3 in WT1). For details of the strains, see Materials 
and Methods. The data collapse (Fig. [5j\, orange circles) corresponds well to the predicted curve (dashed 
red line). Second, the activity-dependence of the demethylation rate is diminished according to Eq. [B] 
when considering adaptation without feedback through activity-dependent CheB phosphorylation, while 
keeping the steady-state activity constant (Fig. [5p, green line). To validate this prediction, a mutant 
strain was created, which contains non-regulatable CheB with about 10 percent of CheB-P activity. 
The CheB expression level was increased to produce the kinase activity of WT2 (A* ml/ 2). All other 
chemotaxis proteins are expressed as in WT2 cells. We find that the experimental rate of FRET-activity 
change from time-course data (green squares) is consistent with this prediction. 

The statistical significance for each of the two predictions (Fig. [5j\ and C) was tested as follows: For 
each prediction, we randomly permuted a number of data points from the control experiment and the 
prediction-testing experiment. Then we calculated the distribution of squared errors between the rates 
of activity change from the model and FRET measurement for the permuted data sets (Fig. [5]B and D). 
For four permuted pairs of data points the error is always above the error for the unpermuted data sets 
(Fig. [5]). For fewer permutations the error lies at the lower bound of the distribution (not shown). This 
confirms that the control and prediction-testing data sets arc significantly different and match our model. 

Discussion 

Sensing and adaptation are fundamental biological processes, enabling cells to respond and adjust to their 
external environment. Adaptation extends the range of stimuli a sensory pathway can respond to, while 
its dynamics determines how long a stimulus will affect the cell's behavior. In this work, we developed a 
model to quantitatively describe experimental dose-response curves, as well as time courses of chemotaxis 
signaling in adapting wild- type cells. Our model includes (i) the signaling activity of two-state mixed 
chemoreceptor complexes in response to added or removed attractant concentration step changes based on 
the Monod-Wyman-Changeux model, (ii) the kinetics of the ligand concentration in the flow chamber, and 
(m) a detailed mechanism for adaptation, including multiple layers of feedback regulation and imprecise 
adaptation. In particular, we find that the finite ligand flow speed and fast, activated demethylation 
explains for the first time the gradually reduced amplitudes in removal dose-response curves for increasing 
ambient concentrations (Fig.[l]). Our adaptation model introduces a strong receptor-activity dependence 
of the demethylation rate, and hence is able to reproduce the observed asymmetry of slow adaptation to 
addition of attractant and fast adaptation to removal of attractant (Fig.|2p). Such dynamics yields long 
runs up the gradient and short tumbles, sufficient for random reorientation of the cell and escape from 
potential toxins. Furthermore, this strong activity dependence has the additional benefit of reducing the 
fluctuations in receptor methylation level introduced by the adaptation mechanism itself. We found for 
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the total variance of the receptor-complex methylation level (<5M 2 }=0.87 compared to 2 for a previous 
model for precise adaptation with weaker activity dependence (details of the calculation can be found in 
the Supplementary Text SI). This is because a fluctuation in the receptor methylation level leads to an 
increased change in activity and hence increased rate to return to the adapted activity. 

Our model for precise adaptation predicts the data collapse of adaptation time-courses, allowing the 
convenient study of the adaptation dynamics (Fig. EK). Specifically, the data collapse allows to evaluate 
the effects of ligand flow and adaptation dynamics, as well as imprecise adaptation. We found that 
adaptation to large concentration step changes is significantly imprecise (sec Supplementary Text SI for 
further details). We also extracted the kinetics of the receptor methylation level dm/dt from experimental 
time courses from the data collapse (Fig. [3)3) , which is difficult to measure directly when relying on the 
quantification of the receptor methylation level using standard biochemical methods [251136) . According 
to our model, the activity-dependence of the receptor methylation level is a monotonously decreasing 
function of the receptor complex activity. Ultimately, this kinetics determines the compromise between 
long memory of previous concentrations and quick recovery for sensing new concentration changes |14) . 
Furthermore, we experimentally tested two predictions to validate our adaptation model. We analyzed 
the effect on the adaptation dynamics when changing the adapted receptor activity, as well as introducing 
a non-regulatablc CheB mutant to remove the negative feedback from phosphorylation of ChcB by the 
kinase ChcA. In both cases, our model is consistent with experimental measurements (Fig. [5]), supporting 
the finding of multiple layers of feedback regulation in adaptation. 

While the MWC model is relatively well established [2"4ll4"lT(4"4"] , we also considered alternative models 
for receptor signaling. These include a phase-separation model with mixed complexes separating into ho- 
mogeneous complexes of Tar and Tsr at high ambient concentrations, as well as a lattice model with finite 
coupling between neighboring receptors (see Supplementary Text SI). Lattice models were previously 
suggested 54,55 , including a lattice formed by coupled CheA molecules [56], but were found to be incon- 
sistent with FRET data [ST]- We found that the dynamic MWC model describes dose- response curves 
far better than the alternative receptor signaling models investigated, particularly the reduced response 
amplitudes upon removal of attractant. Furthermore, the data collapse we introduced in this paper en- 
abled us to compare different adaptation models proposed in the literature with FRET time-course data 
(Fig. H]). We found that while several models are consistent with the data, our model compared most 
favorably with the data. 

We chose a simple model for adaptation with very few fitting parameters to explain the observed 
asymmetry in adaptation time-courses, i.e. slow adaptation to addition and fast adaptation to removal 
of attractant. Compared to the static MWC model, there are minor discrepancies between our model 
and experimental addition dose-response curves (Fig. [T]). However, these can be rectified by refitting 
the dynamic MWC model by adjusting adaptation rates and receptor complex size simultaneously (see 
Supplementary Text SI), or by choosing an adaptation model with a more complex activity dependence. 
It should also be noted that adaptation rates needed to accurately describe dose-response curves are 
larger than those found when fitting the adaptation dynamics to the data collapse. This discrepancy may 
in part be due to using only a single set of experimental data for the data collapse, while dose-response 
curves were averaged over at least three sets. In addition, more complex processes not taken into account 
in our simple adaptation model, e.g. limited supply of substrate (methionine) for methylation, or the 
binding and unbinding kinetics of ligand, may be important for describing the dynamics. 

Although our adaptation model is independent of biochemical details, our predicted fast demethy- 
lation at high activities may be due to cooperativity of two ChcB-P molecules. According to in vitro 
experiments, CheB-P binding to the carboxyl-terminus of chemoreceptors induces an allosteric activation 
of the receptor, increasing the demethylation rate [J0J. However, in contrast to CheR, CheB-P binds 
only weakly to the tether [3D]. Reconciling these two observations, it is conceivable that two CheB-P 
molecules are necessary for efficient demethylation at high activities: one CheB-P molecule may bind to 
a tether to allosterically activate a group of receptors (assistance neighborhood), while another CheB-P 
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molecule demethylates the receptors in the neighborhood. As two ChcB-P molecules are not required to 
bind to the same receptor, this mechanism is not contradicted by the small number of ChcB molecules in 
a cell. An alternative, simpler mechanism for cooperativity is dimerization of ChcB-P molecules, which, 
however, has not been observed experimentally [221158] . 

Our adaptation model likely also applies to attractants other than McAsp, since the dynamics of 
adaptation only depend on the activity of receptor complexes, independent of the details of external 
ligand concentration. According to the MWC model, different attractants (or their mixture) are inte- 
grated at the level of the free-energy of a receptor complex, which determines its activity. However, the 
imprecision of adaptation we found in FRET time courses at large MeAsp concentrations is in contrast to 
earlier experiments, which showed that the frequency of tumbling adapts precisely to aspartate, but not 
serine |28U29j . The imprecision in adaptation to serine is readily explained if the number of Tsr receptors 
is larger than the number of Tar receptors per complex, since the available receptor methylation sites 
in a complex are more quickly used up in response to serine binding to Tsr receptors [42U49| . However, 
the ratio of Tar and Tsr per complex is strongly dependent on the growth conditions, making a defini- 
tive conclusion difficult [5 9) . Future experiments may show if the imprecision observed in adaptation to 
McAsp in FRET is reflected also in the tumbling frequency, or if imprecise adaptation is compensated 
for in order to yield perfect adaptation at the behavioral level. 

Materials and Methods 

Strains Two different wild-type strains of E. coli were used. Wild-type strain 1 (WT1) is VS104 
A(cheY cheZ) that expresses the FRET pair consisting of CheY-YFP (YFP; yellow fluorescent protein) 
and its phosphatase CheZ-CFP (CFP; cyan fluorescent protein) from a pTrc-based plasmid pVS88, which 
carries pBR replication origin and ampicillin resistance and is inducible by isopropyl /3-D-thiogalactoside 
(IPTG) [23]. Wild-type strain 2 (WT2) is VS124 A(cheB cheY cheZ) transformed with pVS88 and 
pVS91, which carries pACYC replication origin and chloramphenicol resistance and encodes wild-type 
CheB under control of pBAD promoter inducible by L-arabinose. The ChcB-mutant strain is VS124 
A(cheB cheY cheZ) transformed with pVS88 and pVS97, which is identical to pVS91 except it encodes 
the non-reg ulatable CheB D56E . The D56E mutation was introduced into CheB by PCR. It prevents ChcB 
phosphorylation, but allows protein to retain basal level of activity. Cells were grown at 275 rpm in a 
rotary shaker to mid-exponential phase (ODgoo ~ 0.48) in tryptone broth (TB) medium supplemented 
with 100/ig/ml ampicillin, 34/ig/ml chloramphenicol, and 50 fxM IPTG. WT and ChcB mutant strains 
were induced by and 0.0003% arabinose, respectively, to produce comparable kinase activity (as assessed 
by the change in the level of FRET upon saturating stimulation). The ChcB protein level was estimated 
using Western blots with CheB antibodies, and was at approximately 0.5-fold (WT2) and approximately 
5-fold (CheB D56E mutant) the native level of CheB. 

FRET measurements The experimental procedures follow those detailed by Sourjik and Berg |23j . 
Cells were tethered to a cover slip, and placed in a flow chamber. Cells were subject to a constant fluid 
flow of buffer or MeAsp at indicated concentration (flow speeds 1000 /d/min for WT1, and 500 fA/mrn for 
WT2 and CheB mutant, respectively). Concentration step changes were achieved by abruptly switching 
between buffer and MeAsp, or different McAsp concentrations. Fluorescence resonance energy transfer 
(FRET) between excited donor, CheZ-CFP, and acceptor, phosphorylated CheY-YFP, in a population 
of 300-500 cells was monitored using an epifluorescence microscopy setup. Emission light from CFP 
and YFP was collected and their intensity ratio R was used to calculate the time-dependent number 
of interacting FRET pairs of CheZ-CFP and phosphorylated CheY-YFP in the cell population, which 
reflects the intracellular kinase activity |23j . The number of FRET pairs normalized by its adapted prc- 
stimulus value (after adaptation to the ambient concentration, but before concentration step changes) was 
calculated from the ratio R according to (R ~ R ) / [{AY / AC) - R] / {{R p[C - R ) / [(AY / AC) - R p[e }} [23]. 
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The parameters R$ and R pre are the ratio for a saturating dose of attractant and the pre-stimulus value, 
respectively, both of which are measured in each experiment. The fluorescence efficiency ratio AY/ AC 
is determined by the experimental setup [50], and was 0.43 (AY/AC=2.3) for WT1 (WT2 and ChcB 
mutant) experiments. FRET measurements were taken with a time resolution of 0.2 s (1 s) for WT1 
(WT2 and CheB mutant). 



Static MWC model This model describes the response of adapted mixed receptor complexes to 
instantaneous MeAsp concentration step changes [24]|43||44] . According to this model, the activity of a 
mixed receptor complex is given by A = [1 + exp(-F)] , where 



F = N 



, l + c/Kf\ , fl + c/Kf 
e(m) + ^ a ln - ' a + v s In ' 



l + c/K° n J ' \l + c/K 



(5) 



is the free-energy difference between the on and off states of the complex. The indexes a and s denote Tar 
and Tsr receptor, respectively. We assumed fractions of Tar and Tsr in a complex according to their wild- 
type ratio, v a :v s =l:lA. The ligand dissociation constants for MeAsp are K° n =0.5mM, K° s =0.02 mM, 
K° n =10 6 mM, and K° =100 mM [13] • The free-energy contribution e(m) is attributed to methylation, 
and was recently extracted from dose-response curves for mutants resembling fixed methylation states [H] . 
We used the interpolating function e(m) = l— 0.5m (for data and fit see Inset of Fig. [2]A_). All energies 
are measured in units of fc^T (fee being the Boltzmann constant and T the absolute temperature). 
Exponential rate constants for the ligand flow were obtained from fits to ligand flow profiles (cf. Fig.[23): 
with A a dd=0.6s _1 and A ram =0.5s _1 for flow speed 1000/d/min, and A a dd=0.27s _1 and A rom =0.28s _1 for 
flow speed 500 /xl/min. The receptor complex size N was estimated from least-squares fits to individual 
addition dose-response curves corresponding to specific ambient concentrations (and therefore adapted 
methylation levels). Note that complex size for removal may be different for each data point as cells are 
adapted to the increased concentration after each step change. The complex size grows with ambient 
concentration [4TU52"] in a roughly linear fashion, N(cQ)—a a +aiCQ with a =17.5 and ax=3.35 mM . Both, 
individually fitted N values, as well as the fitting function N(c ), are shown in Fig. [21V. We assumed 
an adapted receptor complex activity A* = 1/2.9 ~ 0.34 for WT1 as assessed from the maximum and 
minimum values of the experimental dose-response data in Fig. [TJ Steady-state activities for WT2 and 
ChcB mutant were estimated to be A* ml/2. For comparison of model and data, we normalized the 
receptor-complex activity for WT1, WT2 and ChcB mutant by their respective activities when adapted 
to ambient concentration. 



Precise adaptation The dynamic MWC model combines the static MWC model with a dynamical 
model for adaptation. In our model for precise adaptation, the rate of change of the average receptor 
methylation level m is given by (Eq. [5]) 

dm , .. ,o 
— =g R (l-A)-g B A 3 . 

The methylation and demethylation rates do not depend directly on the concentration of MeAsp or the 
methylation level, only on the receptor complex activity A. Such dynamics leads to precise adaptation 
of the receptor complex activity to adapted level A* for a constant MeAsp stimulus j42]|45]. This 
model takes into account the receptor-activity dependence of the methylation and demethylation rates, 
as well as additional layers of feedback regulation for demethylation by CheB, including the activation 
of demethylation enzyme CheB by phosphorylation and potential cooperativity between phosphorylated 
ChcB molecules. For Fig. [T][3j we determined the demethylation rate <7b=0.11 s" 1 from a least-squares 
fit to addition and removal dose-response curves (WT1) using the receptor complex size N(cq) from the 
static MWC model. The methylation rate ^=0.0069 s _1 is given by the condition that at steady state 
(dm/dt=0) the activity equals A*. The fit to the data collapse in Fig. 2] resulted in (?^=0.0019 s _1 (and 
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5b=0.030 s -1 ), used in Fig. @] and [5] for WT1. For WT2 in Fig. [SJA., we used the same methylation rate 
constant as for WT1, but adjusted the dcmcthylation rate constant to account for the increased adapted 
activity A*. For the CheB mutant in Fig. [5p, the rate of change of the average receptor methylation 
level m is predicted to be 

-^ = g R (l - A) - g B A, (6) 

where we assume that the methylation rate is the same as for wild-type cells. The demethylation rate 
constant 5b = <?s^4* 2 =3b/4 includes the basal activity of non-phosphorylatable CheB. Hence, the only 
dependence of the demethylation rate on receptor complex activity is due to binding of CheB to active 
receptors. 



Imprecise adaptation We incorporate the effect of imprecise adaptation, as suggested by time courses 
(cf. Fig. HP), by making methylation and demethylation rates for wild-type cells (WT1) depend on the 
methylation level [49] 

dm m roax -m m 3 

~TT^ gR —^{l-A)-g B —F? A a . (7) 

at m max — m + K m + A 

The parameter m max is the maximum number of methylation sites per receptor, K is the lower bound 
for the number of sites, which need to be available for efficient methylation or dcmcthylation. We use 
m max =4.1 to only allow Tar (not Tsr) receptors to become methylated (the total number of methylation 
sites of a receptor homodimcr being 8). Further, we use K = 0.5 to implement reduced efficiency of 
methylation or demethylation at a low number of available sites. Figure [2p shows time courses for 
adaptation to two concentration step changes using the precise and imprecise adaptation model {gn and 
<7b are the same in both models). The imprecise adaptation model fits the time courses shown far better. 
However, there is a large variability of imprecision seen in different data sets and more experiments are 
needed to produce a general model of imprecise adaptation. 



Rate of activity change To calculate the rate of activity change, the time courses for adaptation 
to step concentration changes were smoothed by averaging every 20 subsequent data points starting 
approximately 10 s after the step onset. The derivative dA/dt was approximated by the difference 
quotient. 
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Figure 1. Response of wild-type cells to step changes Ac of MeAsp concentration at 
different ambient concentrations. Dosc-rcsponse curves: Symbols represent averaged response from 
FRET data (WT1) after adaptation to ambient concentrations 0, 0.1, 0.5 and 2 mM as measured by 
Sourjik and Berg [23] (filled and open circles correspond to response to addition and removal of 
attractant, respectively). Solid lines represent the dynamic MWC model of mixed Tar/Tsr-receptor 
complexes including ligand rise (addition) and fall (removal), as well as adaptation (receptor 
methylation) dynamics. (Inset) Dose-response curves for MWC model without adaptation dynamics 
(lines). FRET and receptor complex activities were normalized by adapted pre-stimulus values at each 
ambient concentration. Squared errors between model and experimental data for the four dose-response 
curves shown are 0.64 (dynamic MWC model) and 3.95 (static MWC model), respectively. For 
comparison, fitting to eight addition and removal dose-response curves using K^fs)' ^afa)> as weu as a 
linear function N(co) as fitting parameters, yields squared errors 0.83 (dynamic MWC model) and 2.09 
(static MWC model), see Supplementary Text SI. 
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Figure 2. Model ingredients. (A) Size of adapted receptor complex N (total number of Tar and Tsr 

receptors per complex) as function of ambient concentration Co (corresponding to adapted methylatioii 
level m). Individual complex sizes (symbols) were obtained by fitting MWC model to dose-response 
curves for addition of MeAsp. These values were fitted by a linear function (line). (A Inset) Energy 
contribution to receptor complex free energy from methylation level m per receptor dimer. Shown are 
fitting parameters for Tar receptors from |41j (symbols), as well as linear fit e(m)=l — 0.5m (in units of 
fcgT with kg the Boltzmann constant and T absolute temperature). (B) Profile of concentration step 
change as measured experimentally using fluorescent marker (solid black line) [23] . exponential fit used 
in dynamic MWC model for WT1 MeAsp profile (gray line; rate constants A a dd~0-6/s and A rem =0.5/s), 
and perfect step change (black dashed line). Addition and removal times are marked by arrows. (B 
Inset) Response of mixed receptor complex to MeAsp removal for perfect (black dashed line) and 
exponentially fitted step change (gray line). (C) Typical time courses of receptor complex activity in 
response to two different MeAsp concentration step changes, Ac=0.05mM (top) and Ac=0.4mM 
(bottom), at ambient concentration co=0.1mM. Experimental FRET measurement (thin black line), as 
well as dynamic MWC model for precise (gray lines) and imprecise adaptation (black lines; m max =A.\ 
and K=0.5). FRET and receptor complex activities were normalized by adapted pre-stimulus values 
before addition of MeAsp. 
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Figure 3. Adaptation dynamics as function of receptor activity for WT1 at ambient 
concentration co=0.1mM for addition and subsequent removal of small (red lines and 
symbols) and large (blue lines and symbols) MeAsp concentration step changes, as well as 
removal of MeAsp back to zero ambient concentration (buffer; green lines and symbols). 

(A) Rate of change of receptor complex activity dA/dt as occurs during adaptation. Thick gray line is 
the analytical result from the dynamic MWC model, where activity change is purely from adaptation 
(methylation) dA/dt=(dA/dm)(dm/dt)=f(A). Colored lines show results from simulated time series for 
small (Ac=0.03 mM) and large (Ac=0.4mM) concentration step changes in MeAsp concentration, with 
activity dynamics recorded starting 10 s after the onset of concentration step change. Precise (solid 
lines), as well as imprecise adaptation (dashed lines; TO ma;E =4.1 and K=0.5) are considered. 
(A Inset) Rate of FRET activity change from experimental time-course data. Small (Ac=0.03 mM) and 
large (Ac=2mM) concentration step changes. (B) Rate of change of the methylation level dm/dt 
corresponding to panel A (normalized by adapted activity A* and C=N/2, where N is the receptor 
complex size) . Effective rate of change of methylation level for all time courses is obtained by 
(dA/dt) /[A(l — A)). (B Inset) Effective rate of change of methylation level from experimental 
time-course data. FRET and receptor complex activities, as well as activity rate changes were 
normalized by adapted pre-stimulus activities at ambient concentrations before addition of MeAsp. 
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Figure 4. Comparison of different adaptation models. (A) Rate of activity change during 
adaptation as a function of activity for FRET data (WT1; symbols) and different adaptation models 
(colored lines). Experimental FRET activity change is measured at ambient concentration co=0TmM 
for added and subsequently removed concentration step changes Ac=0.03, 0.05, 0.1, 0.4 and 2 mM. For 
the five models, the dependencies of the methylation and demethylation rates on the receptor complex 
activity A are given in the legend and are explained in the text. Models were fitted to the experimental 
dA/dt data using a least-squares fit, where the methylation rate constant gn was the only fitting 
parameter. The demethylation rate qb was determined to produce the adapted activity A*«l/3. The 
parameters K\ and K 2 were converted from [17] . (B) Representative time courses for the different 
models in panel A for a concentration step change Ac=0.1 mM at ambient concentration co=0.1 mM. 
FRET and receptor complex activities, as well as activity rate changes were normalized by adapted 
pre-stimulus activities at ambient concentrations before addition of McAsp. Least-squares errors 
between experimental data and model in panel A are 0.0021 [1 - A, A 3 }, 0.0022 [1 - A, A 2 ], 0.0025 
[1 - A, A], 0.0025 [(1 - A) /{I -A + K x ), A 2 /{A + K 2 )], and 0.0036 [(1 - A) /{I -A + K x ), A/ (A + K 2 )]. 
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Figure 5. Effects of (A) steady-state activity and (C) CheB regulation by phosphorylation. 

(A) Black and orange dots correspond to the rate of FRET activity change from experimental 
time-course data for WT1 (Fig. |4j) and for WT2 (addition and subsequent removal of concentration step 
change Ac=0.03mM at zero ambient concentration), respectively. Red lines correspond to the predicted 
rate of activity change dA/dt=f(A) purely from adaptation (solid and dashed lines correspond to 
steady-state activities A* ml/3 and 1/2, respectively). The methylation rate constant ^=0.0019 s _1 is 
the same in each case. Dotted lines indicate bins used to quantify the difference between data sets in 
panel B. (B) Distribution of squared errors (x 2 ) between predicted rate of activity change and 
experimental data sets for WT1 and WT2, when randomly permuting 10 4 pairs of data points between 
the data sets, one pair chosen within each bin in panel A. The error is calculated as the sum of errors 
for each data set (including the permuted data points) against its respective model. The error of the 
unpcrmutcd data sets is indicated by the arrow. (C) Green squares represent the rate of FRET activity 
change from experimental time-course data for CheB mutant (addition and subsequent removal of 
concentration step changes Ac=0.03mM and 0.1 mM at zero ambient concentration). The green line 
represents the rate of change of receptor complex activity purely from adaptation. Orange dots and red 
dashed line are the same as in panel A. Dotted lines indicate bins used to quantify the difference 
between data sets in panel D. (D, left) Distribution of data points of the rate of activity change for 
activities above A=1A WT2 and CheB mutant data in panel C. (D, right) Distribution of squared errors 
between predicted rate of activity change and experimental data sets for WT2 and CheB mutant, when 
randomly permuting 10 4 pairs of data points between the data sets, one pair chosen within each bin in 
panel C. The error is calculated as the sum of errors for each data set (including the permuted data 
points) against its respective model. The error of the unpermuted data sets is indicated by the arrow. 
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1 Review of chemotaxis signaling pathway 

The bacterium Escherichia coli chemotaxes by utilizing a biased random walk towards a nutrient source 
(or away from a toxin source) [IH1]- The swimming path consists of runs, i.e. straight swimming 
driven by coherent motion of flagella, and tumbles characterized by lack of net movement and random 
reorientation of the cell. 

The molecular components of the chemotaxis signaling pathway and relationships between them 
are well-characterized [5], and are shown schematically in Fig. [TJ Transmembrane chemoreceptors 
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Figure 1: Schematics of chemotaxis signaling and its measurement by FRET. (A) Chemotaxis signaling pathway 
in E. coli from receptors to rotary motors and flagella, including phosphotransfcr from CheA to CheY and ChcB, 
ChcY-P diffusion to rotary motors, and dephosphorylation of ChcY-P by phosphatase CheZ. Adaptation involves 
receptor mcthylation by CheR and dcmethylation by CheB-P. (B) FRET pair CFP/YFP used in experiments 
by Sourjik and Berg [B]; YFP tags CheY and CFP tags CheZ. CFP is excited by laser light and transfers its 
energy to YFP during FRET. YFP de-excites by fluorescence. (C) Fast switching between on (white discs) and 
off state (black discs) of an MWC receptor complex. 

localize predominantly at cell poles, where they form large clusters. There are five different types of 
chemoreceptors, each with specific sensing capabilities. The two most abundant receptor types, Tar 
and Tsr, bind respectively the amino acids aspartate (and its non-metabolizable analogue MeAsp) 
and serine. Tsr also binds aspartate and MeAsp with much lower affinity. Binding of ligand induces 
signaling by the receptor across the membrane to the kinase CheA. CheA as well as the protein CheW 
(not shown in Fig. [T]) have been suggested to be involved in receptor-receptor coupling and signal 
integration. When active, CheA autophosphorylates and rapidly passes on a phosphoryl group to its 
response regulators CheY and CheB. Phosphorylated CheY (CheY-P) diffuses to the rotary motors 
which drive the cell's flagella. Upon binding to the motors, CheY-P induces a switch in rotary direction 
resulting in tumbling. CheZ is a phosphatase of CheY-P. Attractant binding reduces the activity of 
CheA, lowering the concentration of CheY-P in the cell, and therefore suppressing tumbling. In 
contrast, repellents cause an increase of activity, enhancing tumbling. 

Adaptation is mediated by the proteins CheR and CheB. CheR methylates receptors to enhance 
their signaling activity. Phosphorylated CheB (CheB-P) demethylates receptors, which reduces their 
activity. During persistent stimulation by a chemical, the combined effect of receptor methylation 
by CheR and demethylation by CheB-P leads to adaptation of the kinase activity to a steady-state, 
allowing the sensing of new changes in attractant or repellent concentrations. 

2 Whole-pathway model 

We consider the following reactions shown in Fig. [TJ (1) auto-phosphorylation of CheA and formation 
of CheA-P (concentrations [A p ]) when receptors are active, (2) phosphorylation of CheY and formation 
of CheY-P (3) association of CheY-P and CheZ ([YJ,Z]), leading to the dephosphorylation of 

CheY-P and dissociation into CheY and CheZ, and (4) phosphorylation of CheB and formation of 
CheB-P ([B p ]). 

Assuming the law of mass- action, our model comprises the following set of ordinary differential 
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equations: 

- A ■ k A ([A] tot - [A p ]) - k Y ([Y] tot - [Y p ]) [A p ] - k B ([B] tot - [B p ]) [A p ] (1) 



dt 



ChcA autophosphorylation CheY phosphorylation CheB phosphorylation 



d[Y p . 
dt 



d[Y p Z] 
dt 



d[B p 
dt 



ky ([Y] to t - [Y p ]) [Ap] - h ([Z] tot - [Y p Z]) [Y p ] + k 2 [Y p Z] 



CheY phosphorylation 



h([Z] tot - [Y p Z])[Yp_ 



CheY-P / ChcZ association 

k B ([B]tot - [B P ]) [Ap] 

V v ' 

CheB phosphorylation 



CheY-P/CheZ CheY-P/CheZ 
association dissociation 

(k 2 + h) [Y P Z] 

V v ' 

CheY-P / ChcZ dissociation 
and CheY-P dephosphorylation 



k- B [B 



(2) 



(3) 



(4) 



spontaneous 
dephosphorylation of CheB-P 



where the ki (with i = 1, 2, 3, A, B, — B and Y) are kinetic rate constants for the individual reactions. 
The activity A of a receptor complex in Eq. [1] is determined by the MWC model, given by 
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l + e 1 



with 
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e(m) + v a In 



off 



1 + c/K™ 



+ v s In 



1 + c 



/Kf \ 



l + c/K° n J 



(5) 
(6) 



as described in the main text. In addition, we include the adaptation dynamics by methylation and 
demethylation of receptors, catalyzed by CheR and CheB-P, respectively (cf. Eq. 2 in the main text), 



dm 
~dt 



g R (1- A) - 

methylation by Che-R 

g R {\ -A)- g B A 3 



i 2 A 



3 b [B p \ 

demethylation by CheB-P 



(7) 



(8) 



where g B and g B are effective rate constants, and g B [B p ] 2 ~ g B A 2 as [B p ] is approximately proportional 
to the receptor complex activity (Fig. [2p). This adaptation model is further explained in the next 
section. The parameter values we used for the whole-pathway model are listed in Table [TJ 



2.1 Rescaling of parameters 

In order to reduce the number of parameters, we normalize the protein concentrations by their 
respective total concentrations in the cell, [A p ] — > [a p ] = [A p ]/[A] tot , [Y p ] — » [y p ] = [Yp]/[Y] tot , 
[Y p Z] — > [y p z] = \Y p Z]/\Y] tot and [B p ] — > [b p ] = [B p ]/[B] tot . Furthermore, we rescale the time by 
the autophosphorylation rate of CheA, k A , t —t t = k A • t, and introduce rescaled rate constants 
according to k\ ->• m = h[Y] to t/kA, k 2 —■ k 2 = k 2 /kA, k 3 -> k 3 = fo/kA, ky -> Ky = fey[^4]tot/^, 
k B — y k b = A;_B[^4]tot/^A and k_ B — > k_ b = k^ B /kA- Overall, this transformation yields dimensionless 
kinetic variables and parameters by measuring phosphorylated protein fractions in units of total protein 
concentrations and rate constants relative to the autophosphorylation rate constant of CheA. Using 
the ratios of total protein concentrations, a\ = [Y] tot /[A] tot , a 2 = [B] tot /[A] tot , and a 3 = [Z] tot /[Y] tot , 
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Figure 2: Steady-state concentrations of individual proteins and CheY-P/CheZ pairs for the whole-pathway 
model for Eq. IMTZl as a function of the receptor complex activity A. Note the different scales of the vertical 
axes. The adapted activity is A* «l/3. 



we obtain the transformed set of equations 

(9) 
(10) 

(11) 
(12) 

The transformed equation for the methylation level of receptors is obtained by replacing time t — > r, 
9R^1R = gnl^A and g B -> 1b = 9BB^ ot /k A in Eq. [3 yielding 





= A-(l- 


[a p ]) - aiK Y (1 - [y p ]) [a p ] - a 2 n B (1 - [b p ]) [a p 




= «y (1 - 


[y p ]) [ap] - Ki (a 3 - [y p ^]) [y p ] + /€ 2 [y p z] 




= Kl(0!3- 


[y P ^])[yp] - ( K 2 + «*) [y P z] 




= «B (1 ~ 


IM) [°p] ~ k -s[^p1- 



dm 
d7 



7i? (l-A)-7i?[&p] 2 A 



(13) 



The new parameter values of this transformed model are listed in Table [TJ 
2.2 Steady-state concentrations 

We analyzed the steady-state concentrations of phosphorylated proteins and CheY-P/CheZ pairs. 
Setting the time-derivatives of Eq. [9lll2lto zero, we solved for the steady-state concentrations of CheA- 
P, CheY-P and CheB-P, as well as the concentration of CheY-P/CheZ pairs as a function of the 
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Table 1: Parameters of the whole-pathway model for chemotaxis signaling for Eq. [T][71 including references 
to literature values where possible, and rescaled parameters for Eq. [91 1131 The literature values are given in 
parentheses where different from our parameter values, ki was determined by the condition that at steady-state 
with A* = l/3, the concentration [Y p ]* = [F]tot/3 [8]. gn was determined by the steady-state activity A* and the 
value for <jb- 



Parameter 


Value 


Reference 


Scaled 
parameter 


Value 


[-4] tot 


5 /xM 


El 


Oi\ 


1.94 


[B] tot 


0.28 uM 

r 


t — i 


«2 


0.056 


Mtot 


9.7 /iM 


M 


«3 


0.113 


[Z]tot 


1.1 fM 




Ky 


50 


k A 


10 s" 1 


® 


K B 


7.5 


ky 


100 iiM" 1 s" 1 




«1 


4.88 


k B 


15 fiM- 1 s' 1 


m 


^2 


0.05 


h 


5.0 fiM^ 1 s- 1 




^3 


20 


k 2 


0.5 s" 1 


m 




0.135 


ks 


200 s" 1 


(30 s^ 1 ) m 




0.0006 


ks 


1.35 s^ 1 


(0.35 s- 1 ) dDE] 


IB 


0.0246 


9R 


0.006 s" 1 








9b 


3.14 /uM- 2 s- 1 









receptor complex activity A. The results are shown in Fig. [2j CheA-P shows a strong non-linear 
dependence on the activity A, i.e., it is strongly activated at high receptor complex activity. It is also 
notable that only a small fraction of the CheA concentration is phosphorylated at maximal receptor 
activity A = 1, which nicely fits estimates from in vitro measurements [7]. All other phosphorylated 
fractions of protein, as well as the concentration of CheY-P/CheZ pairs are approximately proportional 
to receptor complex activity A. 

2.3 Time courses and steady-state assumption 

We tested if the phosphorylation and CheY-P/CheZ association reactions, Eq. l9lfl~2"l are in quasi-steady 
state compared to the slower methylation and demethylation reactions of receptors, Eq. [T3j For this 
purpose, we increased all rate constants for phosphorylation, dephosphorylation, as well as CheY- 
P/CheZ association and dissociation by one order of magnitude, such that concentrations are forced to 
be in quasi steady-state at each time point. Comparing the results to the time courses with the original 
parameter values shown in Fig. [3l we found only minor deviations (exemplified in insets). Therefore, 
the above mentioned reactions are indeed in quasi-steady state to a good approximation. This, together 
with the approximate linearity of the steady-state concentration of CheY-P/CheZ pairs as function 
of receptor complex activity A, permits us to replace the number of FRET (CheY-P/CheZ) pairs by 
the receptor complex activity (with appropriate proportionality factors) as assumed in the main text. 
Similarly, Eq. 2 in the main text arises by replacing CheB-P concentration in the demethylation rate 
in Eq. [13] by the receptor complex activity A, where the methylation and demethylation rate constants 
are g R = j R k A and g B = -f B k A ([b p ]/A) 2 w ~f B kA, respectively. 
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Figure 3: Time courses for the concentrations of phosphorylated proteins and CheY-P/CheZ pairs according 
to the whole-pathway model Eq. IMT51 using ambient McAsp concentration co=0.1 mM and two different con- 
centration step sizes, Ac=0.05 mM (column A) and Ac=0.4 mM (column B). Thick solid curves represent the 
model with parameter values as in table [TJ thin dashed lines assume the quasi steady-state for all phosphory- 
lation reactions (phosphorylation and dcphosphorylation kinetic constants in the model increased by a factor 
20). Insets zoom into the dynamics around the addition or removal time, respectively. Concentrations were 
normalized to their respective adapted values. 
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Figure 4: Dynamic MWC model. (A) Same as Fig. 1 in the main text (main panel), however showing additional, 
previously unpublished data. Shown are dose-response curves for wild-type cells to step changes of McAsp 
concentration (after adaptation to ambient concentrations 0, 0.03, 0.1, 0.3, 0.5, 1, 2 and 5 mM). Symbols 
represent averaged response from FRET data as measured by Sourjik and Berg [6]. Filled and open circles 
correspond to response to addition and removal of attractant, respectively. Solid lines represent the dynamic 
MWC model of mixed Tar/Tsr- receptor complexes, including ligand rise (addition) and fall (removal), as well as 
adaptation (receptor methylation) dynamics. (B) Best global fit of dynamic MWC model with fitting parameters 
.95=0.127 s - \ AT° ff =0.02 mM, AT° n =0.50 mM, Kf=216 mM, K° n =l0 e mM, as well as a = 22 and a 1= 9.6 
mM -1 for the total receptor complex size N = a + (XiCq. 

3 Additional data and best fit of dynamic MWC model 

In Fig. H] we show additional, previously unpublished dose-response data measured as described in 
the main text (cf. Fig. 1 in the main text). The model in panel A is the dynamic MWC model 
from the main text. Panel B shows the best fit of the dynamic MWC model, where we used the 
demethylation rate constant qb , the coefficients determining the receptor complex size as a linear 
function of ambient concentration, and the ligand dissociation constants of Tar and Tsr as fitting 
parameters. We found that parameters overall stay similar to the previously used parameters; in 
particular the ligand dissociation constants do not change significantly. The main difference is larger 
receptor complex sizes than determined by fitting the static MWC model to individual addition dose- 
response curves. To compensate for the larger complex sizes, the adaptation rates are also slightly 
increased, marking the trade-off between increased activity responses by larger complex sizes and 
reduced activity responses by faster adaptation (controlling for MeAsp concentration dynamics). 

Figure [5] quantifies the difference between measured dose-response curves and the static, as well as 
the dynamic MWC model, respectively, as detailed in the main text (cf. Fig. 1 in main text). We plot 
squared errors for each addition and removal dose-response curve. While the error for the dynamic 
MWC model is slightly larger for addition curves, its error for removal curves is much smaller than 
that for the static MWC model. Hence, the dynamic MWC model is suited better to describe the 
experimental data. 

Figure [6] shows the free-energy change associated with each concentration step change. For increasing 
ambient concentrations, the free-energy changes generally decrease at a fixed concentration step change 
Ac. This is the reason for the reduced response amplitudes in the dynamic MWC model at large MeAsp 
step removals, because adaptation compensates for smaller free-energy changes at increasing ambient 
concentrations. 
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Figure 5: Residual absolute squared errors per addition (left panel) and removal (right panel) dose-response 
curve for the static and dynamic MWC model as shown in Fig. 1 in the main text. Note the different axis scales 
for addition and removal plots. 
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Figure 6: Changes in the free-energy difference 5F = F — F* of a mixed-receptor complex upon concentration 
step changes Ac of MeAsp (lines), where F* is the adapted free-energy difference. The curves correspond to 
ambient concentrations co=0, 0.03, 0.1, 0.3, 0.5, 1, 2 and 5 mM with free-energy differences for experimental 
concentration step changes indicated by symbols (filled circles for addition, open circles for removal of MeAsp). 



4 Parameters for static and dynamic MWC model 

In Tab. El we list all parameters of the static and dynamic MWC model used for Fig. 1-3 in the main 
text and Fig. H]A in the Supplementary Text SI, as well as Fig. UJ3 and[HK in Text SI. 

5 Effect of receptor complex size on data collapse 

We found from fitting the MWC model to dose-response curves from FRET that receptor complex size 
increases with ambient concentration (Fig. 2 A in the main text). Hence, we would like to determine 
how the data collapse depends on this effect. According to Eq. 4 in the main text, the rate of activity 
change is proportional to the receptor complex size N. As we do not have a model which describes how 
receptor complex size changes in time in response to concentration changes, we plot in Fig. [7] the data 
collapse for different N corresponding to the concentrations used in the experiments. This provides 
the envelope in which the data collapse is expected to change with N. We find that the data collapse 
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Table 2: Fitting parameters for the static and dynamic MWC model. Parameters include dissociation constants 
for Tar and Tsr receptors in the on and off states, respectively, K° s , K° n , K° s , K° n [13], the parameters of the 
linear approximation of the dependence of the receptor complex size on ambient concentration, iV(co) = ao+aico, 
as well as methylation and demcthylation constants, ga and gs in Eq. 2 in the main text, respectively. Fitted 
parameters are indicated by crosses. 



Parameter 


Fig. 1-3 (main text), HA (Text SI) 


Fig. SB (Text SI) 


Fig. [S]A (Text SI) 


static 




dynamic 


dynamic (best fit) 


static (best fit) 


K° n [mM] 


0.02 




0.02 




0.02 


X 


0.056 


X 


[mM] 


0.5 




0.5 




0.50 


X 


0.15 


X 


Kf [mM] 


100 




100 




216 


X 


100 




K° n [mM] 


10 6 




10 6 




10 6 


X 


10 6 




a 


17.5 


X 


17.5 




22 


X 


37 


X 


ai [mM" 1 ] 


3.35 


X 


3.35 




9.6 


X 


-0.78 


X 


9R [s- 1 ] 


N/A 




0.0069 




0.0079 




N/A 




9B [s" 1 ] 


N/A 




0.11 


X 


0.127 


X 


N/A 





does not change very much compared to the data collapse for ambient concentration cq, and hence we 
neglected the effect of changing complex size in Fig. 3 in the main text. 

6 Unsuitable receptor signaling models 

We tested four alternative models for receptor signaling in an attempt to find a model, which describes 
the gradually reduced response amplitudes upon MeAsp removals at increasing ambient concentrations 
(cf. Fig. 1 in main text) without relying on adaptation and MeAsp dynamics. Dose-response curves 
for each model are shown in Fig. [HJ We found none of the models produced a satisfying fit to the 
experimental data. 

A Saturation model 

While ligand dissociation constants for the Tar receptor were previously determined from FRET 
data [13], slightly different values may lead to saturation of Tar receptors at smaller concentrations 
of MeAsp and reduced response amplitudes. Figure [8]A. shows a fit of the static MWC model to 
addition as well as removal data. We fitted the parameters of the linear relationship between the 
receptor complex size and ambient concentration Co, as well as the ligand dissociation constants for 
the Tar receptor, K° n and K° s . We find an unsatisfying fit, especially for the response to addition of 
MeAsp. Furthermore, the determined receptor complex size decreases with ambient concentration (see 
Inset). This contradicts experiments which indicate an increasing receptor complex size [14], as well 
as stabilization of polar receptor clusters with increasing receptor methylation level (corresponding to 
increasing ambient concentration) [15] . 

B Imprecise adaptation model 

Figure [HJ3 shows the effect of imprecise adaptation on the response amplitudes. For simplicity, we 
assume a linear decline of the adapted activity A*{c$) with increasing ambient concentration Co, with 
A*{Q) = 1/3 and a 20 percent imprecision at concentration 10 mM (see Inset). We observe that 
imprecise adaptation has only a small effect on the response amplitudes. Furthermore, imprecise 
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Figure 7: Effect of ligand concentration-dependent receptor complex size N(c) on the predicted data collapse 
according to Eq. 4 in the main text. Plotted arc the predicted functions f(A) for N corresponding to ambient 
concentration c = 0.1 mM (thick gray line), zero ambient (buffer; black dotted line), and concentration upon 
small (red line) and large (blue line) concentration changes as in Fig. 3 in the main text. 

adaptation tends to increase response amplitudes at high ambient concentrations due to normalization 
by a decreasing value of A*(cq). 

C Phase-separation model 

In this model, a fraction w of mixed receptor complexes composed of Tar and Tsr receptors form 
homogeneous receptor complexes of only Tar and only Tsr receptors at high ambient concentrations. 
This separation reduces activity amplitudes at concentrations below the ligand dissociation constant 
K° s for Tsr, as complexes of Tsr do not contribute to the response. The total activity from mixed 
and homogeneous receptor complexes is 



The individual activities of mixed, A mixed , and homogeneous receptor complexes of Tar, A Ta,r , and Tsr, 
were calculated according to the static MWC model. Mixed receptor complexes are composed of 
Tar and Tsr with ratio v a : v s =\:\A. Homogeneous receptor complexes of Tar and Tsr, respectively, 
have the same ratio. The resulting dose-response curves for this model are shown in Fig. [8)0, assuming 
the probability of separation w(cq) from the Inset. As the ambient concentration does not correspond 
nicely to the data points with decreasing response, we did not find a well-fitting function w(cq). 
Furthermore, this model predicts a smaller response to MeAsp when cells are pre-adapted to a ligand 
for which Tsr, but not Tar is sensitive (e.g. Serine). This contradicts experiments, which show that 
cells remain sensitive (Ref. [16] and V.S., manuscript in preparation). 

D Receptor lattice model 

In the static MWC model the absolute cooperativity of the receptors in a complex results in a saturating 
response upon removal of attractant (cf. Inset of Fig. 1 in main text). Here, we consider an Ising lattice 
of Nt two-state receptor trimers, where each trimer is coupled to neighboring trimers with finite 
interaction strength. We ask if a weaker coupling between receptors can describe the dose-response 
data, and in particular the reduced response amplitudes for removals. 

Figure [8p shows dose-response curves for different interaction strengths. We find, that in order to 
describe the addition data, a strong interaction between neighboring trimers has to be assumed. In 
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Figure 8: Alternative models for receptor signaling. (A) Saturation model. MWC model with optimized 
dissociation constants K° s — 0.056 mM and K° n = 0.15 mM. (A Inset) Total receptor complex size N = ao+aico 
for fitted parameters ao=37 and ai=-0.78 mM -1 . (B) Imprecise adaptation model. (B Inset) The steady-state 
activity decreases with ambient concentration, A*(cq)/A*(Q) = 1 — (0.2/10 mM)co, where 20% imprecision is 
reached at 10 mM. (C) Phase-separation model. Receptor complexes are found in separated complexes with 
probability w depending on the ambient concentration cq. Receptor complex size is assumed constant, ./V =18, 
for mixed and homogeneous receptor complexes. (C Inset) The probability w(co) = pi + P2 Cg c + P3 1 with p 1 =0.01, 
P2=0-99, ^3=0.8 mM. (D) Receptor lattice model. Mixed trimers of Tar and Tsr dimers are arranged on a 
4x4 square lattice with periodic boundary conditions. The average activity of the lattice was calculated by 
exact enumeration. An attractive interaction between neighboring trimers in the same state was assumed, with 
interaction energy J=-0.4 fc^T (solid line), J=-0.35 fc^T (dotted line), and J=-0.3 fc^T (dashed line). (D Inset) 
Corresponding distributions of activities from all states (lattice configurations) when adapted to zero ambient 
concentration. 
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this limit, the lattice model resembles the MWC model where all lattice sites are infinitely strongly 
coupled with all other receptors. In the Inset of Fig. [Hp we show the activity distribution from all 
lattice states. As expected, the distribution becomes increasingly bimodal around the two states with 
all receptors on and all receptors off. 

In the following we describe the details of the model and our simulations. We used a 4-by-4 square 
lattice of mixed receptor trimers with periodic boundary conditions. Each trimer consisted of Tar 
and Tsr receptors with probabilities v a and u s , respectively, where v a :v s =\:\A is the in vivo ratio of 
Tar and Tsr in a cell. The distribution of Tar and Tsr in trimers on the lattice was the same in all 
simulations. Furthermore, each trimer has only two states, on and off. We numerate all possible states 
of the whole lattice (in total n = 2 16 states for a 4-by-4 lattice, i.e. A^=16 receptor trimers). Assuming 
the lattice is in equilibrium, we can calculate the distribution of individual lattice states, and hence 
the average activity of the lattice. The probability of each lattice state depends on its energy, which 
has a contribution from the free-energy difference between the on and off states of each trimer and 
from the interaction between neighboring trimers. The free-energy difference of trimer j is computed 
according to the MWC model 



where the index I = a,s describes the receptor type, Tar or Tsr, within a trimer. The average 
methylation level of receptors in a trimer j is denoted by m? . The methylation energy is e(m J ) = 
3 • (1 - 0.5m j ). 

The interaction energy between neighboring trimers depends on their respective states. If they are 
in the same state (both on or both off), we assign the interaction energy J, if they are in different 
states, we assign the interaction energy — J. The total energy of a lattice state k is determined 
by summing over all free-energy differences of individual trimers and interaction energies between 
neighboring trimers. 

The methylation level m? of each trimer j cannot be calculated analytically due to the finite coupling 
strength between receptor trimers, and hence was determined numerically using our adaptation model 

d ^=g R {l-Ai)-g B {Aif. (16) 
According to this model, the methylation level m? of the trimer j depends on its average activity 



n 



Z 
k=l 

where Z = e~ Ek is the partition function, i.e. the sum over all lattice states and si is the state 
(l=on, 0=off) of trimer j in lattice state k. The steady-state of Eq. 1161 determines the methylation 
level of each trimer, and therefore the adapted free-energy difference e{m 3 ) in Eq. PT5l The average 
activity of the whole lattice is determined by calculating the average trimer activity 

N T 



N T . 



where Nt is the number of trimers on the lattice. 
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7 Comparison of different adaptation models 



In Fig. 4 in the main text we compare different adaptation models to FRET data by collapsing the 
time courses, plotting the rate of activity change dA/dt as a function of the activity A. Here, we 
describe in detail the different models analyzed. In all of the models we assume precise adaptation, 
i.e. that methylation and demethylation rates only depend on the receptor complex activity. For 
each adaptation model, we use a least-squares fit to the FRET data to determine the methylation 
and demethylation rate constants, assuming an adapted activity A* =1/3 and receptor complex size 
N = 17.8. The parameters and quality of fit x 2 f° r each of the models are listed in Tab. El 
Our model Eq. 2 in the main text 

^ = g R (l -A)- g B A* (19) 

is denoted by "(1 — A), A 3 " , referring to the activity dependence of the methylation and demethylation 
rates, respectively. The best fit to the rate of activity change from FRET (fitting parameter <7r=0.0019 
s" 1 , resulting in <7#=0.030 s _1 and quality of fit x 2 =0.0021), and a representative time course for this 
model are shown in Fig. 4A and B in the main text, respectively (red solid lines). Note that this 
model describes the experimental data well, even at high activities. This model also shows a strong 
asymmetry in the time course with slow adaptation to addition and rapid adaptation to removal of 
MeAsp (cf. Fig. 2C in the main text). 

We considered a variation of this model, denoted by "(1 — A), A 2 " , without cooperativity of CheB-P 
molecules, 

^ = g R (l - A) - g B A\ (20) 

where only one CheB-P molecule is necessary for demethylation of a receptor. Together with one 
factor A from the activity of receptors, this leads to a demethylation rate proportional to A 2 . While 
this model is almost as well-suited to describe the rate of activity change from FRET as our main 
model (fitting parameter g/j=0.0031 s _1 ; #b=0.017 s _1 , x 2= 0-0022; see Fig. 4A in main text), the 
asymmetry of adaptation to addition and removal of MeAsp is less pronounced (Fig. 4B in main text). 
Fitting dose-response data using this adaptation model resulted in adaptation rates which were much 
higher than observed in FRET time courses. 

Furthermore, the model denoted by "(1 - A), A" without CheB-P feedback [T9H22] 

d ^ = g R {l-A)-g B A (21) 

yields the fitting parameter ##=0.0048 s , resulting in <7#=0.0091 s _1 and quality of fit x 2= 0.0025. 
Both, the fit of this model to the rate of activity change from FRET, and time courses, are described 
worse than with the other two models. 

Another class of adaptation models was proposed in Ref. [T7] where the idea of ultrasensitivity to 
the adaptation dynamics of CheR and CheB-P was introduced. This model was proposed mainly for 
small changes in activity, such as for fluctuations around the steady-state activity due to noise. We 
denote by "(1 — A)/(l — A + K\), A/(A + K 2 )" the following model 

dm 1 — A A 

-a =9r i-a + k 1 ~ 9b aTk;- (22) 

In this model, CheR (CheB) methylates (demethylates) inactive (active) receptors with Michaelis- 
Menten-type kinetics with Michaelis-Menten constant K\ (K2). In this model, there is no CheB-P 
feedback on the demethylation rate. If K\ and K 2 are small, the adaptation rate depends only weakly 
on the receptor activity. This results in long adaptation (relaxation) times, as well as strong sensitivity 
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Table 3: Parameters of the adaptation models when fitted to the rate of activity change from FRET shown 
in Fig. 4A in the main text. The size of receptor complexes was assumed to be N = 17.8 in all models. a 
Ki = K r /[T] and K 2 = K b /[T], where K r =0.39 [iM and K b =0.54 pM are taken from Ref. [17]. b K 2 = K b /[T] 
with K b =l.25 yM [18]. The concentration of receptors is [T]=17 fjM. 



Adaptation model 


g R (fitted) 


other parameters 


x 2 


"(1 - A), A iv 


0.0019 s- 1 


c/ B =0.030 s- 1 


0.0021 


U {l-A), A 2 " 


0.0031 s" 1 


5B =0.017 s" 1 


0.0022 


a 1— A A n 
l-A+Ki ' A+K 2 


0.0188 s" 1 


95=0.020 s" 1 
Ki= 0.0229 a 
K 2 = 0.0318 a 


0.0036 


a 1-A A 1 » 
l-A+Kx ' A+K 2 


0.0046 s- 1 


95=0.014 s" 1 
Ki= 0.0229 a 
K 2 = 0.0318 a 


0.0025 




0.00318 s- 1 


95=0.014 s^ 1 
K 2 = 0.0735 6 


0.0032 



to protein fluctuations of either CheR or CheB through rates gR and 95. We used K\ = K r /[T] = 
0.0229 and K 2 = Kh/[T] = 0.0318, where we took K r and Kb from [17] and the concentration of 
receptors is [T]=17 (jM.. As shown in Fig. 4A in the main text, the model without CheB-P feedback 
"(1 -A)/(l — A + K t ), A/ (A + K 2 ) n does not describe the rate of activity change from FRET (fitting 
parameter 95=0.0188 s ; 95=0.020 s -1 , x 2 =0.0036). Furthermore, the time course shown in panel B 
looks qualitatively different from experimental time courses (cf. Fig. 2C in main text). 

A variant of the model also includes CheB-P feedback, which introduces another factor A in the 
demethylation rate 53 . We denote this model by "(1 — A)/{1 — A + K\), A 2 /(A + K 2 )" , which 
corresponds to 

dm I- A A 2 

-dF =9R l-A + K l - 9B ATK- 2 (23) 

This model fits the FRET activity change in Fig. 4A in the main text relatively well (fitting parameter 
95=0.0046 s -1 ; 95=0.014 s" 1 , x 2 =0-0025). However, this model is not very different from the simpler 
model "(1 — A), A 2 " , as the CheB-P feedback introduces a strong activity-dependence. 

In the model suggested by Barkai and Leibler [18] CheR methylation does not depend on the activity 
state of receptors, and hence active, as well as inactive receptors get methylated. The kinetics of the 
methylation level is described by 

dm A 

-dT = 9R - 9B ATK- 2 > (24) 

where the parameter value K 2 = i^/[T] =0.074 with Kb=1.25 /jM. [18], and [T] as above. Note 
that this model is a special case of above model "(1 — A)/(l — A + Ki), A/ (A + K 2 )" with Ki=0. 
Fitting to the FRET activity change yields 95=0.00318 s , resulting in 95=0.014 s _1 and quality 
of fit x 2 =0.0032. The predicted data collapse, as well as time courses are very similar to the model 
"(1 - A) /(I - A + Ki), A/ (A + K 2 y\ and is therefore not plotted in Fig. 4 in the main text. 

8 Analysis of adaptation noise 

The receptor methylation level is subject to fluctuations due to the random nature of methylation and 
demethylation events. However, the adaptation dynamics also filters fluctuations in ligand concentra- 
tion (translated into fluctuations of the receptor activity) , averaging over and smoothing high-frequency 
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noise by its slower dynamics. Here, we estimate the variance of the methylation level of a receptor 
complex due to these two noise sources. Equation 2 of the main text describes the deterministic kinetics 
of the average methylation level of receptors in a mixed receptor complex, 

^ = g R (l-A)-g B A 3 . (25) 

Now, we consider the kinetics of the total methylation level of a receptor complex. The total methy- 
lation level M is the sum of the individual methylation levels rrii of all receptors in a complex, 
M = Y^i=i m ii with N the number of receptors per complex. The rate of change of the total methy- 
lation level is 

^ = N R k R (l-A)-N B k B A 3 , (26) 
at 

where we explicitly indicated the number of the modifying CheR and CheB-P molecules, N R and Nb, 
respectively. The modification rates for a single receptor are related to the rates for a receptor complex 
via g R = N R k R /N and g B = N B k B /N, respectively. To describe fluctuations about the mean total 
methylation level due to methylation and demethylation events, we introduce the noise rj(t) and write 

^i = N B k R (l-A)-N B k B A 3 + v (t). (27) 
at 

We assume r](t) is the sum of individual noise terms contributed from each modifying enzyme CheR 
and CheB-P acting on groups of receptors, so-called assistance neighborhoods 



N R Nb 
i=l i=X 

where rj R u\ and i]B(i) are independent Gaussian white noises with zero mean (r] R ^(t)) = (??B(i)(i))=0, 
autocorrelations (vR(i)(t)VR{i)(^)) = QR ' _ and ( r lB(i)('t)''lB(i)('t')) = QB • S(t - t'), and vanishing 
cross-correlations. To estimate the noise intensities q B and q&, we assume that the number of methyl 
groups, which are added (removed) by each enzyme molecule CheR (CheB-P) in a time interval, are 
Poisson distributed, i.e. their variance equals the mean number of added (removed) methyl groups. 
Therefore, the noise intensity q B associated with each CheR molecule is determined by its mean rate 
of methylation, 

q R = k R (l-A*). (29) 
Similarly, the noise intensity qB for demethylation is 

q B = k B A* 3 , (30) 

where we only consider noise from one molecule of CheB-P. We are interested in the steady-state 
fluctuations of the total methylation level. Therefore, we linearize Eq. [26] around the steady state to 
obtain the kinetics of the deviation 5M from the mean methylation level 

( N R k R + 3N B k B A* 2 ) 5 A + n(t) (31) 



dt 

N R k R + W B k B A* 2 ) (j^) ■ SM + ^ • + V (t). (32) 

In the second step, we used that the receptor complex activity is subject to fluctuations from the 
methylation level, as well as the ligand concentration. The derivative of receptor complex activity with 
respect to the free-energy difference (at steady state) is given by 

r)A 

_ = _ A *(1-A*). (33) 
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The total methylation level of a receptor complex enters the free-energy difference through 

2 a V 1 + c/ K T J \ 1 + C I K T J 

=Ef =1 (i-|-i) 

where mi are the methylation levels of receptors i. Therefore, the derivative of the free-energy difference 
F with respect to M is given by 

OF 1 



dM 2' 

The derivative of the free-energy difference F with respect to c is given by 



(35) 



<9c a V c + -^a ff c + K° n J s \c + Kf c + K™ 
In summary, the kinetics of <5M is determined by 

d ( 6M ) _ / jy-^ + 3 NB k B A* 2 ) A*(l - A*) ■ ( i<5M - //<5c ) + r/(t). (37) 



tft V a n J v ' V2 

V v ' V 

=,\ 

To calculate the variance of the methylation level, we Fourier-transform Eq. [37] 

ioo5M = -A Q<5M - /x5c^ + f), (38) 

where the hat symbol denotes the Fourier transform. The power spectrum Sm of fluctuations in M is 
defined as the average of the absolute value squared of 5M 

W-^- S^S . (39) 

Here, denotes the noise intensity of methylation and demethylation, and X 2 /j?(\5c\ 2 ) is due to the 
uncertainty from the ligand concentration^, where we assumed the two contributions are independent. 
In this formula, we see explicitly the noise filtering of fluctuations in ligand concentration by the 
kinetics of the methylation level, given by the frequency-dependent factor. 

In the following, we calculate the variance of the methylation level of a receptor complex only due 
to methylation and demethylation events. As rj(t) is composed of independent white noises, its total 
noise intensity qu is the sum of the individual noise intensities, 

QM = (\fl\ 2 ) = N R q R + N B q B = 2N R q R = 2N R k R (l - A*). (41) 

The last equality uses the fact that at steady state methylation and demethylation rates balance each 
other in Eq. [26] To calculate the variance of the methylation level we need to integrate the power 
spectrum over all frequencies co, 

/AA/f 2\ f <2M 2q M 

< 5M ^ = J 2^^TAV4 = ~T' (42) 



1 Fluctuations of the ligand concentration characterized by (5c 2 ) can be quantified as presented in Ref. [24U25| by 

(5c 2 ) = -5L_ . C) (40) 
naDr 

which corresponds to the time-averaged low- frequency limit of the noise power spectrum |25I26| . The parameter a is the 
size of the ligand binding site of a receptor, D is the ligand diffusion constant, and r is an averaging time due to slower 
downstream reactions. The parameter a is of the order one and depends on further receptor details [251126] . Using 
a ~ 1, a=l nm, D=100 /im 2 /s, a typical ligand concentration c = \J K° a K° n = 0.1 mM [5T], and r = I/Zca = 0.1 s 
corresponding to slow autophosphorylation of CheA, we obtain (5c 2 ) = 5 ■ 10~ 6 mM 2 . 
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and obtain 



= 7 = 7 29R 9N =- A r FT (43) 

{N R k R + W B k B A* 2 )A*(l-A*) {g R + 3g B A* 2 )A* A* + 3(1 -A*) v ' 

= 0.87. 

Here, we used that the adapted activity is A* ~ 1/3, and that the relation between the methylation 
and demethylation rate constants g R and g B is given by the steady state of the methylation kinetics 
Eq.[25j 

I- A* 

A* 

This result can be compared to results for other adaptation models previously reported in the 
literature. Reference [20] uses a linear dependence of methylation and demethylation rates on the 
receptor activity, instead of the nonlinear dependence in Eq. | 



9B = 9R— ^3-. (44) 



^ = g R (l - A) - g B A. (45) 

In an equivalent approach using assistance neighborhoods as described above, the authors calculate 
the variance of the total methylation level to be 

Hence, the variance of the total methylation level of a receptor complex is reduced for adaptation 
kinetics with strong activity dependence of the demethylation rate (Eq. [25]) . compared to the linear 
adaptation model (Eq. |4"5|) . The reason for this is the stronger negative feedback, leading to the rapid 
attenuation of fluctuations in the receptor complex activity. Mathematically, the prefactor of the 
linearized demethylation rate in Eq. [J3] leads to the reduction of the variance of the methylation level 
of the receptor complex. 



9 Quantification of adaptation imprecision 

In Fig. [9] we quantify the imprecision of adaptation. Cells were adapted to 100 /iM ambient con- 
centration with adapted pre-stimulus activity A* re measured by FRET. Concentration step changes 
of various sizes were added, and cells adapted to the new concentration with post-stimulus adapted 
activity ^4p OSt . We define a measure of imprecision as 

A* _ j\* 

Imprecision = — post — . (47) 

We find that adaptation is highly variable from experiment to experiment (high standard deviation). 
However, cells are found to consistently adapt imprecisely at high concentrations. 
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Figure 9: Imprecision of adaptation. FRET time courses were measured for cells adapted to 0.1 mM ambient 
concentration, and subject to various concentration step changes Ac. Levels of adapted FRET activity were 
determined before and after each added concentration step change, and the imprecision was calculated as (^4p OSt — 
A* Te )/A* Ie . Symbols correspond to mean values of imprecision, and error bars indicate the standard mean error 
based on three replicates. The star indicates statistically significant difference from zero with Student's t-tcst p- 
value smaller than 0.05. (Inset) Example FRET time course for Ac=2 mM with adapted pre- and post-stimulus 
activity indicated. 
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